A Primal-Dual Proximal Algorithm for Sparse Template-Based Adaptive Filtering: Application to Seismic Multiple Removal

Unveiling meaningful geophysical information from seismic data requires to deal with both random and structured "noises". As their amplitude may be greater than signals of interest (primaries), additional prior information is especially important in performing efficient signal separation. We address here the problem of multiple reflections, caused by wave-field bouncing between layers. Since only approximate models of these phenomena are available, we propose a flexible framework for time-varying adaptive filtering of seismic signals, using sparse representations, based on inaccurate templates. We recast the joint estimation of adaptive filters and primaries in a new convex variational formulation. This approach allows us to incorporate plausible knowledge about noise statistics, data sparsity and slow filter variation in parsimony-promoting wavelet frames. The designed primal-dual algorithm solves a constrained minimization problem that alleviates standard regularization issues in finding hyperparameters. The approach demonstrates significantly good performance in low signal-to-noise ratio conditions, both for simulated and real field seismic data

Status of the GEO600 gravitational wave detector

The GEO600 laser interferometric gravitational wave detector is approaching the end of its commissioning phase which started in 1995.During a test run in January 2002 the detector was operated for 15 days in a power-recycled michelson configuration. The detector and environmental data which were acquired during this test run were used to test the data analysis code. This paper describes the subsystems of GEO600, the status of the detector by August 2002 and the plans towards the first science run

Euclid in a Taxicab: Sparse Blind Deconvolution with Smoothed l1/l2 Regularization

The l1/l2 ratio regularization function has shown good performance for retrieving sparse signals in a number of recent works, in the context of blind deconvolution. Indeed, it benefits from a scale invariance property much desirable in the blind context. However, the l1/l2 function raises some difficulties when solving the nonconvex and nonsmooth minimization problems resulting from the use of such a penalty term in current restoration methods. In this paper, we propose a new penalty based on a smooth approximation to the l1/l2 function. In addition, we develop a proximal-based algorithm to solve variational problems involving this function and we derive theoretical convergence results. We demonstrate the effectiveness of our method through a comparison with a recent alternating optimization strategy dealing with the exact l1/l2 term, on an application to seismic data blind deconvolution.Comment: 5 page

Passive and active seismic isolation for gravitational radiation detectors and other instruments

Some new passive and active methods for reducing the effects of seismic disturbances on suspended masses are described, with special reference to gravitational radiation detectors in which differential horizontal motions of two or more suspended test masses are monitored. In these methods it is important to be able to determine horizontal seismic accelerations independent of tilts of the ground. Measurement of changes in inclination of the suspension wire of a test mass, relative to a direction defined by a reference arm of long period of oscillation, makes it possible to carry this out over the frequency range of interest for earth-based gravitational radiation detectors. The signal obtained can then be used to compensate for the effects of seismic disturbances on the test mass if necessary. Alternatively the signal corresponding to horizontal acceleration can be used to move the point from which the test mass is suspended in such a way as to reduce the effect of the seismic disturbance and also damp pendulum motions of the suspended test mass. Experimental work with an active anti-seismic system of this type is described

A constrained-based optimization approach for seismic data recovery problems

Random and structured noise both affect seismic data, hiding the reflections of interest (primaries) that carry meaningful geophysical interpretation. When the structured noise is composed of multiple reflections, its adaptive cancellation is obtained through time-varying filtering, compensating inaccuracies in given approximate templates. The under-determined problem can then be formulated as a convex optimization one, providing estimates of both filters and primaries. Within this framework, the criterion to be minimized mainly consists of two parts: a data fidelity term and hard constraints modeling a priori information. This formulation may avoid, or at least facilitate, some parameter determination tasks, usually difficult to perform in inverse problems. Not only classical constraints, such as sparsity, are considered here, but also constraints expressed through hyperplanes, onto which the projection is easy to compute. The latter constraints lead to improved performance by further constraining the space of geophysically sound solutions.Comment: International Conference on Acoustics, Speech and Signal Processing (ICASSP 2014); Special session "Seismic Signal Processing

Fiber-optic push-pull sensor systems

Fiber-optic push-pull sensors are those which exploit the intrinsically differential nature of an interferometer with concommitant benefits in common-mode rejection of undesired effects. Several fiber-optic accelerometer and hydrophone designs are described. Additionally, the recent development at the Naval Postgraduate School of a passive low-cost interferometric signal demodulator permits the development of economical fiber-optic sensor systems

Re-Analysis of HFT Data Using the Apollo Lunar Surface Gravimeter Data

Introduction: The Apollo Passive Seismic Experiment (PSE) was carried out on Apollo 12, 14, 15 and 16. Network observations of four seismic stations were performed for five years from 1972 to 1977. The PSE was a successful mission that informed us of the lunar crustal thickness and seismic velocity structure of the Moon from direct observations of the lunar interior (e.g. [1]). However, the paucity of seismic stations and the limited number of usable seismic events have been a major problem of lunar seismology. An additional observation point enables us to expand the network and the observable area will expand accordingly. Using a data set called the Work Tape, Kawamura et al. (2008) [2] showed that the Lunar Surface Gravimeter (LSG) on Apollo 17 functioned as a seismograph. With this additional seismic station, we tried the first seismic analysis using the LSG data

Time-frequency representation of earthquake accelerograms and inelastic structural response records using the adaptive chirplet decomposition and empirical mode decomposition

In this paper, the adaptive chirplet decomposition combined with the Wigner-Ville transform and the empirical mode decomposition combined with the Hilbert transform are employed to process various non-stationary signals (strong ground motions and structural responses). The efficacy of these two adaptive techniques for capturing the temporal evolution of the frequency content of specific seismic signals is assessed. In this respect, two near-field and two far-field seismic accelerograms are analyzed. Further, a similar analysis is performed for records pertaining to the response of a 20-story steel frame benchmark building excited by one of the four accelerograms scaled by appropriate factors to simulate undamaged and severely damaged conditions for the structure. It is shown that the derived joint time–frequency representations of the response time histories capture quite effectively the influence of non-linearity on the variation of the effective natural frequencies of a structural system during the evolution of a seismic event; in this context, tracing the mean instantaneous frequency of records of critical structural responses is adopted. The study suggests, overall, that the aforementioned techniques are quite viable tools for detecting and monitoring damage to constructed facilities exposed to seismic excitations

Seismic study of stellar convective cores

It has been shown that a discontinuity in the derivatives of the sound speed at the edge of the convective regions inside a star gives rise to a characteristic oscillatory signal in the frequencies of stellar oscillations. This oscillatory signal has been suggested as a means to study the base of the outer convection zone in low mass stars and possibly the outer edge of the convective core in high mass stars. Using stellar models we show that because of a phenomenon similar to aliasing in Fourier transform, it may not be possible to use this signal to detect the convective core. Nevertheless, it may be possible to determine the size of convective cores using the frequency separation \nu_{n+1,l}-\nu_{n,l}.Comment: Accepted for publication in A &